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Abstract 

The renormalization group functions are calculated in D = 4 — e dimensions for 
the ^-theory with two coupling constants associated with an 0(iV)-symmetric 
and a cubic interaction. Divergences are removed by minimal subtraction. The 
critical exponents i], z/, and u are expanded up to order e 5 for the three nontrivial 
fixed points 0(iV)-symmetric, Ising, and cubic. The results suggest the stability 
of the cubic fixed point for N > 3, implying that the critical exponents seen in the 
magnetic transition of three-dimensional cubic crystals are of the cubic universality 
class. This is in contrast to earlier three-loop results which gave N > 3, and thus 
Heisenberg exponents. The numerical differences, however, are less than a percent 
making an experimental distinction of the universality classes very difficult. 
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1 Many experimentally observable features of critical phenomena are described correctly 
by a scalar quantum field theory with a ^-interaction. Field theoretic renormalization 
group techniques [f] in D = 4 — e dimensions [2, 3, 4] combined with Borel resummation 
of the resulting e-expansions [5] have led to accurate determinations of the critical expo- 
nents of all 0(N) universality classes. An extension of the (9(iV)-symmetric theory by 
an interaction E 8 <^ makes the theory applicable to magnetic systems in which the 0(N) 
symmetry is broken and the magnetization is restricted to point along the edges or the 
diagonals of a hypercube in N dimensions. The extended theory interpolates between 
an (9(iV)-symmetric and a cubic system. The (9(iV)-symmetric and the cubic fixed point 
interchange stability depending on N. For N < N C} the first is stable, for N > N C} the 
second is stable (see Fig. f). 



Figure f : The Stability of the fixed points in the ^-theory with (9(iV)-symmetric and 
cubic coupling for N < N c and N > N c . Our results are compatible with N c = 3. 

From previous work based on three-loop calculations it is known that N c must lie 
somewhere between 3 and 4 [6, 7, 8]. Thus, for N < N C} there is symmetry restora- 
tion: Although 0(N) symmetry is broken by the second interaction, it is restored by 
the fluctuations at the critical point. The physically most interesting value of N is 3, 
where the (9(A^)-symmetric fixed point characterizes the critical behavior of the classical 
Heisenberg model of magnetism. The previous result 3 < N c < 4 implied that all mag- 
netic systems with cubic symmetry occurring in nature show an (9(3)-symmetric critical 
behavior. 

The five-loop results to be presented in this letter suggest that the critical value 
N c lies below N = 3, so that the cubic fixed point may govern the critical behavior of 
magnetic transitions in cubic crystals. Due to the vicinity of the Heisenberg fixed point, 
however, the difference in the critical exponents will be hard to measure experimentally. 
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It is well known now to calculate analytically the Feynman graphs of the </> 4 -theory 
in 4 — e dimensions using dimensional regularization [9] and the minimal substraction 
scheme [10] in the five-loop approximation [11]— [17] . The e-expansions of the renor- 
malization group (RG) functions were given in [18] up to e 5 . Semiclassical considera- 
tions have revealed the large-order behavior of the perturbation series [19] and of the 
e-expansions [20]. The two informations together yield the presently known critical 
exponents via several resummation techniques [5]. 

2 In this paper we calculate the five-loop e-expansions for a </> 4 -theory with two coupling 
constants with the euclidean action 



A = I d D x 



^dn<f>idn<f>i + )pn 2 B 4>i4>i + ^-(gi B T} ikl + g 2g T 2 kl ) fafafafa 



where i = 1, 2, . . . , N is a real iV-component field in D = 4 — e dimensions and 
m -B? 9\bi 9zb are the bare mass and coupling constants, respectively. The tensors in 
the two interaction terms have the following symmetrized form: 

T ijki = ^ijhi + SikSji + SuS kj ), (2) 

2 I 1 i i = J = k = h 

T i 3 ki = hjki = \ (3) 
I , else. 

The action is positive definite within the classical stability wedge bounded by the lines 
9i + 92 >0, g\ > and Ng\ + g 2 >0, g\ < 0. This is changed drastically by the fluctu- 
ations (see Fig. 1). 

The renormalization conditions are: 

lf(p) ~ ^Hp)S tJ , (4) 
ri ^ T^ki + ^2 ^ Tf lkU (5) 



r (4) 
1 ijki 



S.P. 



where S.P. denotes the symmetry point of the four-point function at which the incoming 
momenta p 1 ^, p 2 ^, p 3 ^, p 4 ^ satisfy p a p 13 = (4<5"^ — 1). Condition (4) implies that the 
N field components cf>i have only a single wave function renormalization constant. This 
is a consequence of the symmetry of the action under reflection cf>i — > —cf>i and under 
permutations of the N field indices i. The same symmetry guarantees the condition (5) 
to be fullfilled to all orders, 



It should be noted that for N = 2,3, the combinations of T} jkl and Tf ]kl exhaust 
all tensors of rank 4 for which the theory has only one length scale [21]. For N > 4, 
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more tensors are admissible without introducing new length scales [22]; but we shall 
only consider the above two tensors, for simplicity. 

There are four renormalization constants Za (A = </>,2,4i,4 2 ) which relate the bare 
mass mg and two coupling constants §i B of Eq.(l) to the corresponding physical param- 
eters by 



m B = -TT™ 2 = Z m 2m 2 ; g lB = // Z J\ 2 g l for i = 1,2, (6) 



where fi is a mass parameter introduced for dimensional regularization. The RG- 
functions are introduced in the usual way: 

Pi(9i,92) = ^df,g t \ giB , g2B , mB ,e = ^d^g l , (7) 

B 



(8) 



7m{9i,92) = ^d^\ogm\ giB)g2B , mB) e = //d^logm . (9) 

B 

The calculation of the contribution of each graph to the RG-functions proceeds as de- 
scribed in [18] and references therein. For each graph, we evaluate the corresponding 
contractions of the tensors T} jkl and T? w , decomposing the result again with respect to 
these tensors. 

After some work, we obtain the following expressions for the RG-functions up to five 
loops: 

M9i,92) =-eai+ai 2 ^+3132 2+ gi 3 [-N-f]- gi 2 g2 f -31 32 2 1+ 

n 4 l"/V 2 11 , Af / 461 1 20 <(3) \ 370 , 88 <(3) ] 3 r 79 N , 659 , 64 ((3) ] 2 2 [ N , 107 ,o ^/o\l , - - 3 7 , 

31 [W j^+N { + — + — + — 9 J+31 S2 [-36- + 7g- + — ^J+Si 32 [^ + — +8C(3)j+3 1 3 2 7 + 

. 5 \ N 3_B , at2 ( 395 14 <(3) 10 <(4) 80 <(5) \ AT / 1QQ57 1528 <(3) 124 <(4) 2200 <(5) \ 

L 3888 V 243 9 27 81 )' \ 486 81 27 81 ) 

24581 4664 <(3) 352 <(4) 2480 <(5) 1 
486 81 ' 27 27 J ~ l ~ 

. 4. riv2/_7 C(3)\ | a 7 / 1319 184 <(3) 38 <(4) 400 <(5) \ 15967 4856 ((3) 340 <(4) 2560 <(5) 1 

31 32 [J> ^gj 9 JTJV 9 9 27 ) 81 27 9 9 J ~ l ~ 

3i 3 32 2 [A r ( ^ 35 9 (3) ) 13 5 4 4 33 1456 9 c(3) I 64 ^ 4) 2000 9 <( 5 ) ] 1 

31 2 32 3 [iV(-f + ^)- 4 l I -50C(3)-8C(4)-i^l]+3i 32 4 [- ^ -3 C(3)-6 C(4)] + 

a 6 [jV4f^3 «21) +N 3 ( 6289 , 26 <(3) 2< 2 (3) 7 <(4) 305 <(5) 25 <(6) \ 

aI I V 62208 432 J~ \ 31104 ~ 81 27 24 T 243 81 l~ 

N 2 ( 50531 , 8455 <(3) 59< 2 (3) 347 <(4) 7466 <(5) 1775 <(6) 686 <(7) \ 

I 3888 ' 486 81 54 243 162 ' 27 J "'" 



AT / 103849 1 69035 <(3) , 446< 2 (3) 2383 <(4) , 66986 <(5) 7825 <(6) ir>Ar> ^,^ 

ly \y 972 486 81 54 243 81 1 1 ' ) 

17158 , 27382 <(3) 1088 < 2 (3) 8 
81 ' 81 "T" 27 

„ 5 LfS /^161__17C(3)_C(41X ^2 ( 5 

31 32 I -iv ^ 10368 648 36 jTJv 1 j 



17158 27382 <(3) 1088 < 2 (3) 880 <(4) 55028 <(5) 6200 <(6) 25774 <(7) 
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537437 | 116759 <(3) 3148 < 2 (3) 10177 <(4) 75236 <(5) 24050 <(6) 11564 <(7) , 

~T 07 97 ~T 07 Q7 T 1 M" 



„ 4. 2 I" AT2 ( 1921 , 763 <(3) 1 7 <(4) 5 <(5) \ AT ( 270749 , 9230 <(3) 232 < 2 (3) 4841 <(4) 2045 <(5) 1450 <(6) 245 <( 7) \ 
91 92 U \ 10368 648 36 9 )' \ 2592 81 27 108 9 27 3 / 



1314497 171533 <(3) 1384 < 2 (3) 23105 <(4) 96794 <(5) 25400 <(6) 14210 <(7) . 

+ 27 54 + 27 27 + ' I + 

At 



,3,3 at ( 30277 , 344 <(3) 25 <(4) 208 <(5) \ 2281727 , 37789 <(3) 544 < z (3) 337 <(4) 17444 <(5) 1600 <(6) ooro 
91 92 \ 1296 + 27 ~ 6 + 9 ) + 1296 + 27 ~ 9 ~ 3 + 9 ~ 9 +Z ^ (' > J + 

siW ^(fm.ll^+l^il.ii^ + M 

91 32 5 [iffff + ^|^-18C 2 (3) + 2Mll + 2C(5) + 75C(6)] , (10) 

^2) = -eS2+3 32 2 +4 3l 32-S2 3 f -fli S2 2 f -SI 2 32 ^™+S2 4 + 12 C(3)] + 
W [ i f i +48C(3)]+ 3l 2 32 2 [i|f + 3|5+64C(3)]+ 31 3 32 [-JV 2 ^ +JV (§ + ^) + + 2^] + 
32 5 [- 3 ^-78C(3) + 18C(4)-120C(5)]+ 3l g 2 * [- Ifi -387 C(3)+96 C(4)-600 C(5)] + 
ai 2 S2 3 [w(-i|-^+4C(4))-^fi-724C(3) + 184C(4)-^^] + 

PI 3 3 2 2 [iV 2 (| + ^)+iV(-ff -2^ + 12 C (4)-l^)-ifM_Hl^ + M0|(ll_ 960C(5) ] + 

(-1^+^)+^ (-ife-^+^H 

jv ( - mi - + 6*411 _ ) _ 12S29 _ i 68 c (3) +40 c (4) - + 

92 6 [^+ 7 - 3 ^+45(^3)- 1 -i^+987((5)- 6 -^ + 1323((7)] + 

91 92 b [M|I|07 + 35083 2 C(3) +288c2(3) _ 3697^(4) +5920C(5) _ 2100C(6) + 7938C(7) j + 

SiW [Ar(-^ + 221^_ 17C 2 (3) _l^ + ?^_H^i) 

+ 36|||7I + 125459 <(3) + 2266 < 2 (3) _ 2263 C ( 4 ) + 14328 C ( 5 ) _ "575 <(6) +19404 C(7) ] + 

3 i 3 P2 3 [A^^-^-^ + ii^+A^if^ 

+ 93^7 + 224804 <(3) + 3032 < 2 (3) _ 73018 «4) + 155692 <(5) _ 6300 C (6) + 2391 2 C (7)] + 

a i a 2 \ N 3(^21 91 C(3) , C(4) \ N 2 ( 43295 |tr(3 x 4< 2 (3) 121 <(4) 364 <(5) 50 <(6) \ 

91 92 MV ^20736 1296 "T lg J-fiV I 3HQ4 ^ ^ ' 3 24 ^ 27 9 / 



I 54 r" 81 _ 27 _ 54 + 27 _ 54 + ~ 

1279979 , 784621 <(3) 18154 < 2 (3) 83837 <(4) 275510 <(5) 98975 <(6) 43120 <(7) 
162 + 27 54 + 27 27 +3 



[jV4 (__61 ili(£i + C(liUjV 3 f 3557 151 <(3) 4 <(4) 8 <(5) \ 

L V 15552 972 ~ 108 _/ ~ " V 46656 972 27 ~ 81 ) ~ 

N 2 ( 111217 , 2785 <(3) 92< 2 (3) 1055 <(4) 530 <(5) 950 <(6) 98<(7) \ 
I 23328 ' 243 81 162 ' 27 81 9 J 

»r / 95588 , 15742 <(3) 92< 2 (3) 6592 <(4) 34460 <(5) 14750 C(6) ^ ( ^A 
JV ^ 729 81 81 81 81 81 T^»U ) J 

, 389095 , 259358 <(3) 13288 < 2 (3) 27166 <(4) 179696 <(5) 63500 <(6) 28420 <(7) 1 

729 ' 243 ' 81 81 81 81 9 J ' V / 



72(51,52) = 3 l 2 ^+ 3 l 32 i+ 3 2 2 ^- 3 l 3 [|g + |Tl+^]- 3 l 2 3 2 [f + i 
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„ 4 [ 5 N 3 , 5 N 2 , 85 N , 125 1 ,„ 3 „ T 5 N 2 , 5 iV , 125 1 . 2 „ 2 [5 iV , 145 1 . „ 3 65 ,„ 4 65, 

01 5 r r i3-+<L?i) +N 3 ( i^_m) +N 2 ( _M59_ + iMm_isL±i) + 

yL L V 62208 ^ 3888 ^ V 93312 972 J~ \ 11664 ^ 972 162 ) ~ 



4 „, \n 3 ( 65 I 5 C(3) "I I N 2 ( 185 5 <(3) )lN( 395 , 20 CO) 25<(4) \ 1505 , 115 CO) 55 <(4) 1 
yi y2 [JV ^ 20736 T 1296 J " V 7776 216 J " V 216 T 81 54 J 243 T 162 27 J T 

3_ 2 [ 325 JV 2 |V / 4453 , 23 CO) 5 <(4) \ 58177 , 191 CO) 130 <(4) 1 
yL y2 L 31104 ~ \ 3888 ~ 216 27 ) 3888 T 108 27 j~ 

91*92* [iY(-^L + iifl)-I|Iil + -^l-5C(4)] +3 i^ + +32 . [-^iiOl.COl] ,(12) 



7m(9l,92) = 91 ^+92 h-91 2 ™-Pl 92 §- 3 2 2 ^ + 

fll 3 [^i + ili, + +gi 2 g2 Sr^N +gi ^2 25^^3 1 + 

01 4 [jyS «ll) +N 2 ( 947 4 CO) 5 <(4) \ AT / 5777 22 CO) 16 CO) \ 7765 34 CO) 22 CO) 1 

yL L V 15552 108 J~ \ 1944 81 54 J~ V 1944 81 27 ) 1944 81 27 J 

a S \ N 2(_k C(3)\ |V / 421 10C(3) 10 CO) A 7765 68 C(3) 44 C(4) 1 

yi y2 [JV V 1296 9 y 72 27 9 J 324 27 9 J T 

2 2 [ V / 1841 CO) C(4)\ 9199 13 CO) 26 C(4) 1 
i/1 y2 [JV V 864 6 3 j 216 3 3 J T 

W [iv(-|| + ^)-^-i^-6C(4)] +3 ^ [-m-M3)_Mil] + 

01 5 U4 f 7 , 17C(3) C(4)\ , N 3 ( 2831 , 487 C(3) C 2 Q), 23C(4) 5 C(5) 25 C(6) ' 
yL L V 124416 ~ 7776 432 J~ \ 23328 ~ 3888 243 T 1296 486 ~ 486 



^2 I 291907 , 1261 C(3) 149C 2 (3) 437 C(4) 2 CO) , 1475 C(6) , 
1 93312 972 486 648 243 972 ' 

N I MMIi 



5832 1 324 

49477 , 1327 CO) 194 C 2 Q) , 667 C(4) 8 CO) , 775 C(6) " 



Sl 4 S2 [N 3 ( 41 3 4 5 72 I 8 2 5 C 9 ( 2 3) 5 1 C 4 ( 4 4) ) I - 2 ( 113135 I 1175 <(3) 5C - 2(3) I 145C(4) 25C( - 5) I 125C(6) 

N ( ^ + 

J Y I ma T 



108 1 6 162 1 

247385 , 6635 CO) 485 C 2 Q) , 3335 C(4) 20 C(5) 3875 CO) " 



3 2 1^2 / 2045 , 785 C(3) 7 C(4) 7 CO) \ , AT ( 362281 , 3743 CO) 61 C 2 (3) 1493 C(4) 1 1 C(5) 775 C(6) \ 
yi yi ^ 20736 ~ 1296 72 54 J~ \ 10368 ~ 324 27 ~ 216 18 ~ 54 / 



a, 2 a, 3 \N f 26173 , 71 CO) 2C 2 Q) 25C(4) 14 C(5) 25 C(6) \ 
i/1 y2 [JV I 2304 T 48 9 T 12 9 ^ 9 J 

, 32003 , 627 CO) 403 C 2 Q) , 475 C(4) 59 CO) , 3325 C(6) 1 
i ~96~'~4 9 i ~6~ l ~9~'~18j~'~ 

01 4 [ W 26171 77 CO) 7C(4) 2 C(5) \ 589141 , 959 CO) 45C 2 Q) 643 C(4) 19 CO) , 375 C(6) 1 
yi W2 |_J> ^ 13824 96 ~ 16 3 ) ~ 3456 T 12 2 T 16 T 6 T 4j T 

5 [ 158849 , 1519 CO) 9C 2 Q) 65 C(4) CO) , 75 C(6) 1 /I Q\ 

92 L 4608 T 96 2 T 8 T 2 T 4j' V J ! 



3 The three critical exponents rj } v and u completely specify the critical behavior of the 
system [4]. They are determined from the RG-functions of the ^-theory by going to the 
infrared- stable fixed points (g^^g^)- These are found from the zeros of the /3-functions 
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in the form offree e-expansions: 

<?* = = >>r**- (14) 



CO 

k=l 



Besides the Gaussian fixed point (g^= g\ = 0), there are the following three nontrivial 
fixed points, Ising {g\ = §^g\ = g^), symmetric or Heisenberg (g\= g\ -,g\ = 0) and cu- 
bic {g\ = gi \g\ = g%). For N = 1, the Heisenberg fixed point coincides with the Ising 
fixed point (gf = g^), whereas cubic (g^ = —g?) and the Gaussian fixed point will have 
degenerated critical exponents. For N = 2, we will find degenerated exponents for the 
cubic and the Ising fixed point (—g% = g^)- As N increases, the cubic fixed point ap- 
proaches the Heisenberg fixed point from below, crossing it at N = N c , where g% changes 
its sign. For N — > oo, the cubic fixed point moves into the Ising fixed point. 

Since the e-expansions for the Ising and Heisenberg fixed points can be deduced from 
the series in Ref. [18], we exhibit the expansions only for the cubic fixed point 

r< ('-106 + 125 iV-19 N 2 ) 

f 3 [ 22472 45080 , 38329 41971 1955 , / 56 28 8 , 8 \ f(r>\ \ i 

L 729 JV 5 729 N l ~ l ~ 972 N 3 5832 N 2 5832 N ^ \ 9 N i 9 N 3 3 N 2 T 9 N ) "> \° > J ^ 



[-i 



2987 



9683 N 7 6561 N 6 6561 N 5 19683 JV 4 52488 N 3 104976 N 2 314928 N 

' 29680 i 46696 12764 7934 , 2744 , 110 \ 

k 243 N 6 243 N 5 243 JV 4 243 N 3 243 N 2 243 N ) ^ v ' 

' 32 16 16 , 2 \ f/A\ j I 80 , 200 , 80 , 80 80 \ f 

,9JV 4 9N 3 9 JV2"I"3JV JU^T^ 3 N 5 + 27 N i + 27 N 3 + 9 N 2 27 N ) ^ 



I" 176746j 
L 531441 



N 9 531441 N 8 531441 N 7 1062882 N 6 8503056 N 5 

5795035 



" 8503056 JV 4 68024448 N 3 136048896 N 2 136048896 N 



^2187JV 8 729 N 7 2187 N 6 1458 N 5 4374 JV 4 4374 N 3 34992 N 2 T 34992 N 
2213 47369 , 3929 , 217 ^ ^^-j 



IC(3) 



243 N 6 486 N 5 81 JV 4 1944 N 3 486 N 2 ~ 648 N , 

26296 340 



-. 27 N 7 729 N 6 729 N 5 27 JV 4 729 JV 3 729 JV 2 243 JV 



IC(5) 



100 \ ^C R \_i_/ —1078 , 1225 2450 , 539 343 , 98 \ 
27WMW"H qjV 6 + 3W 5 9JV 4i_ 9JV3 9 N 2 + 9 N ) *> <■ > 

W~~sTW^~sT1^^~sTn 3 ~sin 2 ~2Tn) ( 2 ( 3 )] + 0[ e ] > (15) 



81 JV5 T 162 JV 4 ^ 81 N 3 ^ 81 JV 2 2 7 iV" y'V /^y g jy6 3 jy"5 9 JV 4 ^ 9 N 3 9JV 2 T 9JV 
3136 112 2834 , 3089 , 626 296 16 



C> ,„ ■> ('424-534 JV+93 JV 2 + 17 N 3 ) 

i)2 3 N 81 N 3 + 

fZ \ 89888 1 187528 123707 , 90281 , 11713 , 709 , / 224 , 16 , 80 32 i_ \ 

I 2187 JV 5 " 1 " 2187 JV 4 2187 N 3 ^ 8748 N 2 ^ 1 7496 N ^ 1 7496 ^ 2 7 JV 4 ^ 3 N 3 ^ 27 N 2 27 N 27 ) U°,!J T 



p 4 I 23820320 69389720 , 25018256 35478331 , 11944655 



59049 N 7 59049 JV 6 19683 N 5 59049 JV 4 1 18098 N 3 157464 N 2 944784 N ~ 944784 



0450 56 106 \ s(n\ 

. 729 N 6 729 N 5 ^ 243 JV 4 T 81 N 3 729 N 2 81 JV 729 J ^> 

136 
. 27 JV 4 



136 1 28 1 64 28 2_ \ j- 1 *\ , / 320 400 440 80 , 280 , 40 \ ^/pr\ 1 

i "9JV3" l "27JV 2 27 N 27 ) ^ \^>^ \ 9 N 5 27 N i 81 N 3 g N 2 + 81 N 81 J >.°^J i " 
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I" 7069870976 , 2937809504 4398801284 , 9923276525 5033294725 
L 1594323 N 9 177147 N 8 1 77147 N 7 531441 N 6 708588 N 5 

4069429 321451 



2834352 JV 4 ^ 17006112 N 3 22674816 N 2 ^ 45349632 N 408146688 

17618048 , 46032704 40407016 , 10668718 , 4840987 1176529 , 64261 , 10361 11221 \ > 

6561 N 8 ^ 6561 N 7 6561 N 6 ^ 6561 N 5 ^ 13122 JV 4 6561 N 3 ^ 26244 N 2 ^ 104976 N 104976 J > 

72941 323 



. 729 N 6 729 JV 5 162 JV 4 ^ 243 JV 3 5832 JV 2 972 JV 5832 , 

202864 10202 83693 , 29770 , 1628 , 373 



i 29770 i 1628 i 373 \ s ( r\ 
' -vna w2 729 /V" 729 / 'V / 



. 81 JV 7 729 N 6 729 JV 5 729 JV 4 729 JV 3 729 N 2 T 729 JV T 729 , 

' 3100 1525 425 1025 , 275 , 25 \ /l f l , / 4312 14896 , 3626 2254 , 980 98 49 \ 

.81JV 5 81JV 4 162JV 3 81 N 2 54 JV 54 ) ^W^W N 6 27 N 5 " + " 9 JV 4 27 N 3 " + " 27 N 2 S« 2 7 ^ V ^ 



4354 

" 81 N 7 81 JV 6 27 N 5 81 JV 4 ' 



8ff3+^2+^ + ^)C 2 (3)]+0[ e 6 ] . (16) 



The stability ol the fixed points is determined by the eigenvalues u\ and cu 2 of the matrix 

dA( gl ,g 2 ), 

If the real parts of both eigenvalues are positive, the corresponding fixed point is infrared 
stable. The Gaussian fixed point is doubly unstable. At the Ising fixed point, one 
eigenvalue u\ is negativ. The Heisenberg and the cubic fixed point interchange stability 
for N = N c , the former being stable for N < N c where < 0. The stability wedges 
of the critical theory are visible in Fig. 1. They differ from the bare stability wedge. 
Outside these wedges, the transition is of first order [23]. 

To find the crucial number N c determining which fixed point governs the critical 
behavior in D = 3 dimensions, we study the eigenvalue uj^' of the stability matrix as 
a function of N (the other eigenvalue lo^ remains positive and can be ignored), whose 
e-expansion reads 



2 3N y ' 



2 (-848 + 660 TV + 72 N 2 - 19 N 3 ) 

e 77^- tt; (18) 



81 N 3 (2 + TV) 



3 n=*C 3 N* , 4 m ClN* , 5 S^CfiV* 



+ <_> t — U t <± t — U t O 



8748 iV 5 (2 + N f 944784 iV 7 (2 + N) 5 102036672 N 9 (2 + N)\ 



+ 0[t 



The coefficients C\ are listed in Table 1. From the vanishing of uj^' we extract the 
e-expansion of N c 

2 / 5 , 5C(3)\ , 3 / 1 , 5C(3) , 15((4) 25((5) 



iV c = 4-2 e + e ^-- + ^J +e ^-- + ^ + ^-^J + (19) 

e 4 /_ ^_ 93((3) _ 229( 2 (3) 15 ((4) _ 3155 ((5) _ 125 ((6) 11515((7) \ 5 
V 384 128 144 32 1728 12 384 J [€ 1 ' 

The same expansion is found from the condition g% = 0. Unfortunately, this expansion 
is badly divergent making it difficult to calculate the value of N c at e = 1. With the 
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help of Pade approximations we obtain the values 
Pade [1,1]: N c = 3.128 Pade [2, 2]: N c = 2.958 

Pade [2,1]: N c = 2.792 Pade [1,2]: N c = 2.893 

Pade [3,1]: N c = 3.068 Pade [1,3]: N c = 2.972 

The approximation Pade [1, 3] is unreliable since the Pade denominator has a pole at pos- 
itive e, where the exact result should be regular. The highest symmetric approximation 
is usually the most accurate one, from which we deduce the estimate: 

N c « 2.958. (20) 

Before our work, the e-expansion for N c was known only up to the order e 2 [7] so 
that only the Pade [1, l]-approximation was available [7] which yielded N c 3.128. 
An alternative method, the scaling field method of [24] in which the RG-equations 
are transformed into an infinite set of differential equations, was used in Ref. [25] to 
estimate N c directly in D = 3 dimensions yielding N c = 3.4. Thus the previous results 
yielded iV c -values larger than 3 implying that the critical behavior of magnetic systems 
with cubic symmetry is governed by the Heisenberg fixed point. The symmetric Pade 
approximation to our expansion suggests that N c lies below three so that the cubic fixed 
point is the relevant one. 

By inserting the expansions (15) and (16) into (12) and (13), we find the critical 
exponents rjc and vq of the cubic fixed point via the defining relations [4] 

r, = 2 72 ((/*,(/*) , 1/u = 2(1 - j m (gl,g$)) . (21) 

The exponents are 

\jv° = 2+ (JV-1) + [ Ii ^3 (424-326 JV+19JV 2 )] + 

£ 3 [ 7 1 n5 (-359552 + 573728 JV-264936JV 2 + 28358JV 3 +937JV 4 + 4 2 ( ^ 4 2) (-14+11 JV-JV 2 ) C(3))] + 
£ 4 [ 18895 1 68 N T (381125120-923268480 JV+798088608JV 2 -284926360JV 3 +32693424JV 4 + 768780JV 5 + 24857JV 6 ) 
+ g N e (ll8720-152032JV+29816JV 2 + 17936JV 3 -4124JV 4 -119JV 5 ) C(3) 

+( f^(- 14+llw - w2 ) c ( 4 )+ir^( 36 - 2JV - 4Jv2 - 8Jv3+Jv4 ) (: ( 5 )] + 

£ 5 [ 1 q('-452471742464+1470211004416JV-1869697955840JV 2 + 1160186503168JV 3 

L 204073344 N 9 V 

-350446218272 JV 4 +41122747144JV 5 + 144762448JV 6 -68383472 JV 7 +64327JV 8 ) 
+ 7^J-p- (-70472192 + 153589248 JV-106996288JV 2 + 18129888JV 3 +6458072 JV 4 -1726592 JV 5 -8716 JV 6 -3335 JV 7 ) C(3) 
+ 1944 1 n6 (ll8720-152032JV+29816JV 2 + 17936JV 3 -4124JV 4 -119JV 5 ) C(4) 
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+ 2187 N? (- 915840 + 89 7560Af-53320Af 2 -2676Af 3 -86879iV 4 +20128iV 5 + 528 N 6 ) C(5) 
+ 5^5 (36-2JV-4JV 2 -8JV 3 +JV 4 ) C(6) + ^g (44-132 JV+65 JV 2 -9 JV 3 +8 JV 4 -JV 5 ) C(7) 
+ ^- 8 ^ T (-784+588JV+232JV 2 -163 JV 3 -23 JV 4 + 5 JV 5 +JV 6 ) C 2 (3)] [ + 0[e 6 ] . 



(22) 




+ i^^( 28 - 6JV - 16Jv2+4 Jv3 " jv4 ) c ( 3 )] + 




-731940728 JV 4 +46139232JV 5 + 1948700JV 6 +321511JV 7 ) 
+ i749 1 6 ^ 7 (-189952 + 266624JV-31584JV 2 -80376JV 3 + 29704 JV 4 -2196 JV 5 -329 JV 6 ) C(3) 

+ iT^( 28 - 6W - 16w2+4w3 - w4 ) c ( 4 )+?i?^(- 36+Jv2+18Jv3 - 5Jv4+Jv5 ) (: ( 5 )] f + °[ e6 ] 



(23) 



As anticipated, the exponents for the cubic fixed point, are degenerate with those of 
the Gaussian fixed point, taking free-field values. At N = 2, the exponents take Ising 
values. 

For a reliable extraction of the critical exponents from our expansions it will be 
necessary to find the large-order behavior of the perturbation expansions in §\ and g 2 . 
A combination of Borel and Pade resummation methods are expected to give the best 
results. Such methods should confirm the result N c < 3 and render the precise differences 
between the critical exponents of the Heisenberg and the cubic fixed point. From the 
above rough Pade resummations, we expect the differences to be less than a percent. 
This will make it difficult distinguishing the two fixed points experimentally. 
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Co 


2876416 


cl 


22268920832-5633556480C(3)-1074954240C(5) 


cl 


-1740544+580608C(3) 


cl 


-5244000000+5418233856C(3) + 188116992C(4)-2120048640C(5) 


cl 


-3188544+829440C(3) 


C\ 


-21313343616+8161855488C(3)+456855552C(4)-2508226560C(5) 


c\ 


2340592-62208C(3) 


cl 


11104506624-370593792C(3) + 295612416C(4)-1425807360C(5) 


c\ 


54656-404352C(3) 


cl 


1087015200-2657484288C(3)-83980800C(4)+432967680C(5) 


cl 


-162696-114048C(3) 


C'l 


-1171729344-359023104C(3)-173000448C(4)+901393920C(5) 


c% 


15700+7776C(3) 


cl 


54071304+211455360C(3)-67184640C(4)+343388160C(5) 


cl 


-937+2592C(3) 


cl 


26339632 + 27133056C(3)-5878656C(4) + 20062080C(5) 


c* 


-12196003840 


°10 


-3320774-2164320C(3) + 1469664C(4)-9797760C(5) 


c\ 


5267640320-4923555840C(3) 


c» 


-24857+154224C(3) + 209952C(4)-933120C(5) 



C\ 57916383035392 

C\ -14984240431104+35071472959488C(3) 

G\ -143277741441024+39325680009216C(3) + 2022633897984C 2 (3) + 10938734346240C(5) 
C\ 1211449212928-58923130945536C(3) + 5995664769024C 2 (3) 

-1595232092160C(4) + 22263974461440C(5)-2085841207296C(7) 
C\ 162125071257600-92450858139648C(3) + 5046265184256C 2 (3)-3420504391680C(4) 

+9499251179520C(5)-580475289600C(6)+3887249522688C(7) 
C\ 38426107425792 + 21006186430464C(3)-2702434959360C 2 (3)-468572553216C(4) 

-11373660831 744C(5)-1725301555200C(6) + 14600888451072C(7) 
C'l -162464030196224+69899709652992 C(3)-6924425232384C 2 (3)+3943630872576C(4) 

-19961756909568C(5)-2644387430400C(6) + 13700184293376C(7) 
Cf 53925927185664+4311047245824C(3)-3110702579712C 2 (3) + 2963245731840C(4) 

-11389462659072 C(5)-2410584883200C(6) + 5736063320064C(7) 
Cf 14422298978304-19789519380480C(3)+993580204032C 2 (3)-319987003392C(4) 

+ 1995616714752C(5)-874744012800C(6) + 592568524800C(7) 
C| -7438725755776-2538882164736C(3) + 1199003959296C 2 (3)-1007366492160C(4) 

+4815696457728C(5) + 528071270400C(6)-1724374407168C(7) 
Cf -147321611712 + 2076145468416C(3) + 251458670592C 2 (3)-263068176384C(4) 

+ 1308067024896C(5) + 730632960000C(6)-2005844456448C(7) 

258021138336+211760262144C(3)-41520107520C 2 (3)+48221775360C(4) 

-269590685184C(5)+317951308800C(6)-960701720832C(7) 
C\ 2 -8428414672-68664506112 C(3)-17716589568C 2 (3) + 25217754624C(4) 

-141247120896C(5) + 51900134400C(6)-170733806208C(7) 
C\ z -2039889852-1386414144C(3)-851565312C 2 (3) + 1546506432C(4) 

-7514135424C(5)-3086294400C(6) + 12221725824C(7) 
C\± -138029874+936160416C(3)-25194240C 2 (3)-125341344C(4) + 1017287424C(5)-1826582400C(6)+6666395904C(7) 
C\ & -64327+12966480C(3)-20155392C 2 (3) + 12492144C(4)-49268736C(5)-125971200C(6)+370355328C(7) 



Table 1: The constants appearing in Eq. (19) for uj^' ■ 
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